Origins of anomalous electronic structures of epitaxial graphene on silicon carbide 
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On the basis of first-principles calculations, we report that a novel interfacial atomic structure 
occurs between graphene and the surface of silicon carbide, destroying the Dirac point of graphene 
and opening a substantial energy gap there. In the calculated atomic structures, a quasi-periodic 
6x6 domain pattern emerges out of a larger commensurate 6a/3 x 6v / 3^R30 o periodic interfacial re- 
construction, resolving a long standing experimental controversy on the periodicity of the interfacial 
superstructures. Our theoretical energy spectrum shows a gap and midgap states at the Dirac point 
of graphene, which are in excellent agreement with the recently- observed anomalous angle-resolved 
photoemission spectra. Beyond solving unexplained issues in epitaxial graphene, our atomistic study 
may provide a way to engineer the energy gaps of graphene on substrates. 

PACS numbers: 73.20.-r,81.05.Uw,68.35-p,71.20-b 



Graphene, a carbon allotrope, is a two-dimensional 
hexagonal network of carbon atoms which is formed by 
making strong triangular a-bonds of the sp 2 hybridized 
orbitals [H, El . The 7r-orbitals orthogonal to the hexag- 
onal plane of graphene are responsible for its character- 
istic electronic properties, i.e., a relativistic dispersion 
relation near the Fermi level described by the massless 
free particle Dirac equation []], 0, 0, 0] • Following the re- 
port of successful fabrication of mechanically exfoliated 
graphene on the insulating Si02 surface [3, |j|, tremen- 
dous efforts have been devoted to measure and exploit 
the novel physical properties of graphene [l], E| . 

On the other hand, it has been known for the last 
three decades that, when a wide band gap semicon- 
ductor silicon carbide (SiC) is heated up to 1300°C, 
the mono crystalline graphite forms on the SiC (0001) 

face |H 0, 0, 0, 0, EE lH El • Now > b y fine tunin S 
the growth parameters, a single layer of graphene can 
be grown successfully on SiC [13, LL4| . These researches 
have stimulated interests in resolving fundamental mate- 
rial properties [5|. [yl as well as applying the techniques to 
nanoelectronics [13, E3] , with the merits of precise con- 
trol of the number of layers of graphene 15, lid ~v\ 



and 



a possible large scale production [l4| . 

Epitaxial graphene has demonstrated different physi- 
cal properties compared to exfoliated graphene, exhibit- 
ing many controversial experimental observations El, El, 
II El, EE E3, El, El, El, EH . For examnle^scanning tun- 
nelling microscopy (STM) images [l2|, [13, LUl, LL8|, LL9[ show 
a 6 x 6 hexagonal superstructure with respect to the sur- 
face unitcell of 4H-SiC(0001) while low-energy electron 
diffraction (LEED) patterns indicate a larger scale recon- 
struction with a 6\/3 x 6V3R30 periodicity 0, 0, 0, 0, 
0, El, E3, El, El, El El, EE • Moreover, the energy spec- 
trum from STM [l9[ and the angle-resolved photoemis- 
sion spectroscopy (ARPES) measurements [21] show the 
energy gap which still defy precise interpretations [22[. 



Since the potential profiles induced by interfacial atoms 
will play a decisive role in the physical prop erties of 
graphene grown on SiC (0001) 0, 0, EJ El, EE , it is in- 
deed required to know the precise atomic geometries and 
the corresponding electronic structures of the system. 
However, in spite of many experimental observations, the 
atomic structures of the graphene (graphite)-SiC(OOOl) 
interfaces have not been uncovered yet except that they 
have a large scale reconstruction with the 6y/3x6y/3R30° 

periodicity @, 0, 0, 0, 0, El, E3, El, El, E3, El, EE • 

In this Letter, we identify atomic and electronic struc- 
tures of epitaxial graphene on 4H-SiC(0001) by large- 
scale first-principles calculations. In the relaxed atomic 
structures of the 6a/3 x 6V3R30 periodicity with a sin- 
gle layer graphene, a quasi-periodic 6x6 domain pattern 
appears, resolving the aforementioned disagreement be- 
tween LEED patterns and the STM images. The ob- 
tained novel domain pattern originates from interplay 
between strong bonding and lattice mismatch at the 
graphene-SiC(OOOl) interface. With inclusion of another 
layer of graphene, the calculated electronic structures 
show a gap opening and midgap states at the Dirac 
point of graphene, originating from sublattice symme- 
try breaking interactions between graphene and the in- 
terfacial superstructure in the system. Simulated STM 
images and simulated ARPES spectra for the obtained 
atomic structures show excellent agreements with sev- 
eral experimental data [H, El, E3, El, El, El, EH- Our 
study resolved the fundamental issues regarding the role 
of interfaces and substrates in altering physical proper- 
ties of graphene, and thus provides a way to control the 
energy gaps of graphene on substrates. 

We study the atomic and electronic structures of 
graphene and interfacial carbons on 4H-SiC(0001) based 
on ab initio pseudopotential density functional meth- 
ods [11] within the local density approximation [25[ which 
are known to describe the structural and electronic prop- 



FIG. 1: (color online) (a) Side and (b) top views of the atomic 
structure of the buffer layer with the 6 y/3 x 6y / 3^30° peri- 
odicity. The 6a/3 x 6y/3R30° supercell is denoted with blue 
lines and the lxl surface unit cell of 4H-SiC(0001) with 
red. The carbon atoms in the buffer layer are denoted with 
black spheres and the silicon and the carbon atoms in the 4H- 
SiC with orange and green spheres, respectively, (c) Bonding 
characteristics of carbon atoms in the buffer layer. Carbon 
atoms with (without) cr-bonding to surface silicon atoms are 
represented with grey (black) dots. The 7r-bonds, represented 
by black lines, form a super-hexagonal pattern. Four times 
the 6^3 x 6v / 3^30° periodicity is drawn in order to display 
the super- hexagons clearly. The blue arrows are the unit vec- 
tors of the 6\^3 x 6y / 3^30° supercell. 




erties of graphite quite well [2fll2Z|. The 4H-SiC(0001) 
substrate is modelled with four alternating silicon and 
carbon atomic layers, and one or two graphene layers 
are placed on top of the SiC substrate. The atoms be- 
longing to the bottom layer of the slab are passivated by 
hydrogen. To incorporate a large number of atoms in 
the system, we expand the wave function with localized 
basis sets [28]. The basis sets and the pseudopotentials 
are thoroughly tested to reproduce the atomic and elec- 
tronic structures of SiC, graphene and the x V3R30 
model for epitaxial graphene studied in previous liter- 
respectivley. Based on the LEED mea- 

[E B S m El EI EI EI EI 0, the 



atures [2f 
surements 

6V3 x 6^3^30° periodicity is imposed to the SiC(OOOl) 
surface, which is equivalent to 13 x 13 times a graphene 
unit cell [Fig 1(a), (b)]. The atomic positions are deter- 
mined by minimizing the total energy until the forces on 
each atom are less than 0.06 eV/A while atoms belong- 
ing to the last two silicon and carbon layers are fixed to 
the bulk atomic structure of 4H-SiC. To simulate ARPES 
spectra, wavefunctions of the 6\/3 x 6\^3R30° supercell 
obtained with a few k-points are Fourier-transformed to 
the surface reciprocal space of graphene unit cell with 
dense k-points and then integrated along the surface nor- 
mal direction with attenuation corresponding to a photon 
mean free path of 5A. The calculated spectra are broad- 
ened by 30 meV in energy. 

The obtained atomic structure, in the case of one layer 
of graphene on the SiC surface, displays a novel pattern 
of covalent bonding between carbon atoms in graphene 
and silicon atoms on the (0001) face [Fig. 1(c)]. We find 
that, due to the interplay of lattice mismatch and strong 
C-Si bonds, the 6y/3 x 6\^3R30° supercell is split into lat- 
tice matched regions, where carbon atoms are covalently 



FIG. 2: (color online) (a) Simulated spectrum for ARPES 
of the buffer layer on top of 4H-SiC(0001) in the graphene 
Brillouin zone, (b) Simulated STM image of the structure. 
Wavefunctions of which energies lie between the Fermi level 
(Ef) and 0.2 eV above Ef are integrated and the image is 
taken in a plane located at 3 A above the buffer layer. Bright 
(dark) regions correspond to a high (low) current in constant- 
height mode for STM. The 6^3 x 6^/3^30° periodic lattice 
vectors are drawn in blue while 6x6 periodic ones in red. 



bonded to surface silicon atoms with V^3 x \/3R30 o peri- 
odicity with respect to the surface unitcell of SiC(0001). 
Outside the lattice matched regions, there are bound- 
aries consisting of lines of carbon atoms (connected or 
disconnected with each other) without covalent bonding 
to surface silicon atoms. Our calculations show that the 
covalent bonding in the regions cannot be sustained over 
three or four units of the reconstruction and that the car- 
bon atoms with (without) covalent a-bonding to silicon 
atoms move toward (away from) the SiC substrate. The 
small regions of covalent bondings follow the lattice sym- 
metry of graphene so that all the lattice matched regions 
are shown to have large hexagonal shapes approximately 
[Fig. 1]. We find that the novel hexagonal pattern ap- 
pearing in the present simulation is quite robust and does 
not depend on the details of calculations. 

Our atomic structure for a single layer graphene on 
the SiC substrate [Fig. 1] is consistent with several ex- 
perimental observations for the initial stage of graphene 
formation @, H !, EI El EI El, EI El, El , in which 
the observed structure is called a carbon nanomesh [T3] 
or a buffer (or dead) layer [H, EH, 0, [27| . We will call 
it a buffer layer hereafter. First of all, our model for the 
buffer layer has the §^/3 x 6V3R30 reconstruction sat- 
isfying all the LEED measurements [I @, 0, @, i, EI El 



13 EI EI El, Hi Second, our simulated spectrum for 
the ARPES shows that a-bands in the buffer layer ap- 
pear clearly while the linear 7r-bands around the Fermi 
energy are absent [Fig. 2(a)], compatible with the re- 
cent ARPES measurement [20|. Instead, there are several 
fiat bands above and below the Fermi energy originating 
from 7r-orbitals of carbon atoms on the super-hexagonal 
boundaries. In our atomic structure [Fig. 1], the covalent 
bonding between graphene and the SiC surface breaks the 
hexagonal network of 7r-orbitals but preserves a-bonds of 
sp 2 hybridization. Hence, the resulting band structure in 
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the graphene Brillouin zone [Fig. 2(a)] shows no relativis- 
tic dispersion relation near the Fermi level of the system. 
Third, our simulated STM image shows an approximate 
6x6 periodicity with respect to the surface unit cell of 
SiC(OOOl) [Fig. 2(b)]. This also bears a striking similar- 
ity to the observed STM images [H, [H, Q El E, Q , 
although the 6a/3 x 6v // 3^30° periodicity is imposed on 
our atomic structure. The bright regions in the simulated 
STM image [Fig. 2(b)] originate from the (broken) chains 
of 7r-orbitals with approximate super-hexagonal shapes in 
large scale and the dark regions correspond to domains of 
carbon atoms having strong cr-bonds to the surface silicon 
atoms. Our atomic model for the buffer layer has a com- 
plete coverage of carbon atoms on the SiC(OOOl) surface 
without any silion atoms on it. The model is compatible 
with not only the aforementioned spectroscopic measure- 
ment [20| but also recent STM studies on the adsorption 
of metallic clusters [30| and fullerenes [31] on the buffer 
layer indicating clean carbon surfaces. 

Next, we consider another graphene layer on top of the 
reconstructed buffer layer on 4H-SiC(0001) surface [Fig. 
3(a)]. The calculated atomic structure shows an almost 
free standing graphene with corrugations following the 
atomic structures underneath it [Fig. 3(a)]. The corru- 
gation height is ±0.17A, with the mean distance of 3.35 
A from the buffer layer. The simulated STM image [Fig. 
3(b)] shows that all hexagonal networks of carbon atoms 
are clearly visible and approximate 6x6 periodic large 
hexagonal shapes are superimposed on it. This matches 
very well with the existing experimental data [3, El, EH , 
although the 6>/3 x 6V3R30 periodicity is imposed on 
our atomic structure. In our calculations, the apparent 
6x6 periodicity originates both from the atomic cor- 
rugations due to the underlying buffer layer and from 
the weak electronic interaction between graphene and the 
buffer layer. Thus, our results resolve the disagreement 
between the 6x6 periodicity from STM measurements 
and the 6V^3 x 6V3R30 from the LEED measurements. 

We find that the simulated ARPES spectrum of 
graphene on top of the buffer layer show the character- 
istic 7r-bands of graphene as well as a-bands [Fig. 3(c)]. 
The crossing point of two linear 7r-bands of graphene 
(called as the Dirac point [H, Q) is located slightly be- 
low the Fermi energy of the system. However, if looked 
closely, there is a gap opening at the Dirac point. We 
shall defer the discussion of the gap later. From Mul- 
likan population analysis [24], graphene is found to be 
electron-doped with a density, n ~ 8.7 x 10 12 / cm 2 . co nsis- 
tent with experimental observations H, EE E3, 0, EH . 
It is also noticeable that the cr-bands of graphene are 
rigidly shifted up in energy compared with those of the 
buffer layer. This shift is observed in the recent ARPES 
measurement [loj and the size of the shift (1.2 eV) is con- 
sistent with their observations [20]. The energy shift of 
the cr-bands of graphene arises from a potential gradient 
due to a polar nature of the SiC surface. 
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FIG. 3: (color online) (a) Side view of the atomic structure of 
graphene on the buffer layer with the 6V3 x 6A/3.R30 peri- 
odicity, (b) Simulated STM image of graphene shown in (a). 
The unit vectors for the 6^3 x Ga/3^30 periodicity are drawn 
in blue arrows and those for 6 x 6 in red. (c) Simulated ARPES 
spectrum for graphene shown in (a). The bottom of cr-bands 
in graphene is shown to be shifted up in energy (pointed by 
red arrow) compared with that in the buffer layer (blue ar- 
row). The simulated STM image are obtained by the same 
method as in Fig. 2. 



The electronic structure near the Dirac point shows a 
gap of 200 meV with the centre of the gap located at 320 
meV below the Fermi energy [Fig. 4(a)]. It also shows 
midgap states inside the gap. The spectrum exhibiting 
the gap and midgap states is very similar with the recent 
ARPES observations [21]. The gap at the Dirac point 
in our calculation originates mainly from the interlayer 
coupling between graphene and the buffer layer which 
breaks a sublattice symmtery in graphene @, [5l|, [H, Q . 
Considering the valley and pseudo-spin symmetries ex- 
hibited in the electronic structure of graphene []], 0], 
there are two possible ways of inducing a gap at Dirac 
points in a single layer of graphene. One is the mixing of 
electronic states with different pseudo-spins in the same 
valley and the other is the mixing of states belong to 
different valleys 0, Hl|, [H, Q . Due to the presence of 
the buffer layer, there exists a substantially different in- 
teraction at atomic sites belonging to each sublattice of 
graphene [Fig. 4(b)], breaking the sublattice symmetry. 
The resulting wavefunctions near the Dirac point reflect 
such a broken symmetry so that the weight of the wave- 
function on one sublattice is predominant over the other 
[Fig. 4(c)]. The intervalley mixing is found to have a 
negligible contribution to the gap. In case that the sub- 
strate is intentionally removed from the system in the 
calculation, the remaining graphene with present corru- 
gations does not show any gap at the Dirac point. The 
midgap states originate from the interlayer coupling (the 
hopping energy of 0.2^0.3 eV between 7r-orbitals [2[) be- 
tween the 7r-states in graphene and the localized 7r-states 
on the boundaries between super-hexagons in the buffer 
layer [Fig. 4(d)]. The corresponding wavefunctions spa- 
tially spread out into graphene (the topmost surface of 
the system) [Fig. 4(d)] and thus can be detected by sur- 
face sensitive measurements such as ARPES. It explains 
the reason why the energy distribution curve shown in 
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FIG. 4: (color online) (a) The magnified view of the ARPES 
spectrum near the Dirac point, (b) Contour plot (contour 
spacing = 0.1 eV) for the potential generated by the buffer 
layer and SiC substrate only, drawn on a plane located at 
3.35A above the buffer layer. The bright (dark) color corre- 
sponds to high (low) potential. The hexagonal network (grey 
lines) for graphene is drawn to guide the eyes. The dark 
and bright dots represent two sublattices having ~140 meV 
averaged potential difference, respectively, (c) The squared 
amplitude of wavefunctions (isosurface of 3.0 x 10 _4 /A 3 ) of 
which energy is located at the upper apex of the energy bands 
at K (shown in (a)). The amplitude of the wavefunction at 
graphene is denoted in red while the one at the buffer layer 
in light blue, (d) The squared amplitude of the wavefuntions 
whose energies are inside the energy gap at Dirac point (the 
averaged midgap state). The isosurface value is 3.0 x 10 _4 /A 3 . 
The red and blue colours for the amplitude follow the same 
scheme in (c). 



the recent ARPES measurement [2l[ has anomalous non- 
vanishing weights inside the gap at the Dirac point. 

In summary, we have demonstrated that the interface 
between the SiC (0001) surface and graphene show the 
novel large scale atomic reconstruction and the funda- 
mental electronic property of epitaxial graphene is al- 
tered due to the interface. The present atomistic study on 
the atomic and electronic structures of epitaxial graphene 
will play a crucial role not only in designing electronic 
circuits [14| but also in explaining many other puzzling 
observations such as the absence of the quantum Hall ef- 
fect and the weak Shubnikov-de Hass oscillations in the 
high mobility sample 14, 32| . 
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